Regular reflection in self-similar potential flow 
and the sonic criterion 

Volker Elling 

Abstract 

Reflection of a shock from a solid wedge is a classical problem in gas 
dynamics. Depending on the parameters either a regular or a irregular 
(Mach-type) reflection results. We construct regular reflection as an exact 
self-similar solution for potential flow. For some upstream Mach numbers 
Mi and isentropic coefficients 7, a solution exists for all wedge angles 
9 allowed by the sonic criterion. This demonstrates that, at least for 
potential flow, weaker criteria are false. 



1 Introduction 



1.1 The reflection problem 



Reflection of an incident shock from a solid wedge is a classical problem of gas 
dynamics. It has been studied extensively by Ernst Mach |MW75[ IKvdGQl] 
and John von Neumann |Neu43| . as well as many other engineers and mathe- 
maticians. 

Most commonly, reflection is studied in steady inviscid compressible flow, for 
example when shocks in a nozzle are reflected from the walls. The reflections 
can be classified roughly into regular and irregular reflections; see |BD92j for 
a more detailed discussion. In either type, an incident shock Q impinges on a 
solid surface _B (see Figured]). In regular reflection (RR), Q reaches a reflection 
point on the surface, continuing as the reflected shock R (see Figure [T] left) . 



In irregular reflections (IRR), incident and reflected shock are connected by a 
more or less complex interaction pattern which in turn connects to the solid 
surface by a third shock, called Mach stem. The most important irregular 
reflections are double, complex and single Mach reflection (DMR, CMR, SMR); 
various additional types have been proposed |Gud62[ [HBOOl IHT02j . Figure □ 
right shows an (oversimplified) version of single Mach reflection. 
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Figure 1: Left: regular reflection, right: single Mach reflection (oversimplified). 
If the reflection point is steady, the sonic criterion corresponds to Mq > 1. 

The reflection problem has several parameters. For polytropic gas it is sufficient 
to consider the iscntropic coefficient 7 as well as Lq and L/, the Mach numbers 
in the Q rcsp. / regions. The incident shock cannot exist unless Lq > 1. Lq 
and Lj < Lq determine the incident shock (not all Lj may admit a matching 
reflected shock). 

In Mach reflection, the Mach stem, reflected and incident shock appear to meet 
in a triple point. In general this is possible only if they arc joined by a contact 
discontinuity (slip line); for some parameter values it is not possible at all. 
In fact for certain values RR is not possible cither. This is called the von 
Neumann paradox; it is perhaps the most famous of the many problems arising 
in reflection. Many ideas have been proposed towards the resolution of the 
paradox (see e.g. |Gud621 ITRMI IGRT99( IHB00( IHT02| ): no single explanation 
has been accepted widely so far. 

However, this article is concerned with a different question: it is natural to 
ask which parameters cause a RR and which yield IRR. Of course both sides 
of Figure [T] are perfectly valid stationary solutions, so the question has to be 
phrased more carefully. For example: 

1. Which of the two is dynamically stable (e.g. asymptotically stable as a 
stationary solution of the time-dependent problem)? 

2. Which of the two is structurally stable under perturbations like down- 
stream nozzles, wall curvature or roughness, interaction with other flow 
patterns, perturbation of the upstream flow to non-constant with curved 
incident shock, viscosity, heat conduction, boundary layers, noise, slow 
relaxation to thermal equilibrum and other kinetic effects, dissociation 
etc. 

It is not clear whether these questions are really any better than the original one 
— perhaps both sides of Figure[T]are stable. If so, then the new questions would 
merely fail in a less obvious way, as stability is harder to check than existence. 
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Figure 2: A given upstream velocity (/) with possible downstream velocities 
(red curve) for steady shocks with varying normal. The U shock is unphysical; 
K is the strong and L the weak shock. Shocks cannot turn velocities by more 
than the critical angle r*. 



But let us assume for the sake of the argument that the vague problem "does 
RR or IRR occur" can be expressed in some way as a precise mathematical 
question that selects exactly one of the two choices. 

Among the criteria that have been proposed (see [BD92[ Section 1.5]), three are 
most important. The first criterion, called detachment criterion^ states that RR 
occurs whenever a reflected shock exists. Clearly RR is not possible without a 
reflected shock, so this is the weakest possible criterion. 

The velocity vj in the / region of Figure |4] forms an angle t with B; the reflected 
shock must turn this velocity by r so that vj^ is parallel to the wall, satisfying 
a slip boundary condition. 

Given the / region data and 7, let the reflected shock be steady and pass through 
the reflection point, but vary its angle. This yields a one-parameter family of 
velocities vr, forming a curve called shock polar (see Figure [2]). For physical 
shocks there is a maximum angle r^,, between downstream and upstream velocity, 
is determined by the upstream state. 

If the angle r between wall and vj region of Figure S] right is bigger than r* , 
no reflected shock exists. If r = r*, there is exactly one reflected shock. For 
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T < T^, however there are two, called weak reflection and strong reflection. We 
encounter another one of the major issues in reflection: which of these two 
should occur? |EL07| have discussed this question in a related problem. 

The flow in the R region can be supersonic or subsonic. If it is supersonic, 
then waves in the R region cannot travel towards the reflection point. If it is 
subsonic, however, they can reach it and interact with it, potentially altering 
the reflection type. This motivates the second criterion, called sonic criterion: 
RR occurs exactly if there is a reflected shock with supersonic R region, i.e. 
Mach number Lfj > 1. 

On the shock polar (Figure [2]), + indicates the point where Mfj = 1; velocities 
right of it are supersonic, left of it subsonic. Hence there is an angle r+ so 
that for T < T+ the weak reflection L has L_r > 1. For r > t+ however it 
has Lr < 1. The strong reflection K is always subsonic in the R region — so 
the sonic criterion has a pleasant property: only the weak reflection is allowed, 
solving the uniqueness problem. Moreover since t+ < r*, the sonic criterion is 
stronger than the detachment criterion. 

The third criterion is motivated by studying what happens when the parameters 
Lj, Lq are varied so that a transition from RR to IRR occurs. One might suspect 
that the pressure in the reflection point in the R, S regions is continuous and 
does not jump during transition. Then the pressure behind the reflected shock 
in RR and the pressure behind the Mach stem in IRR, a shock approximately 
straight and perpendicular to the wall, must be equal at transition. There is a 
very limited set of L/, Lq, 7 for which this happens; the von Neumann criterion 
(sometimes called mechanical equilibrum criterion) states that the transition 
can occur only at those parameters. 

The von Neumann criterion has various problems. Most importantly, for weak 
incident shocks the pressure behind the Mach stem never matches the pres- 
sure below the reflected shock, so RR should occur in all cases, contradicting 
observations. 



1.2 Self-similar reflection 

Reflection can also be studied in self-similar (sometimes called quasi-steady 
or pseudo- steady) flow. In fact this is advantageous: for finding stationary 
solutions, choosing boundary conditions that yield well-posedness, in particular 
uniqueness, can be rather subtle, as evident from the awkward phrasing of 
the RR-or-IRR question above. For initial-value problems, on the other hand, 
uniqueness is expectecG — or at least a necessary property of any interesting 

^ IE1106||E1105I raise doubt about the Cauchy problem for the Euler equations, but at least 
for potential flow the author expects uniqueness to hold. 
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Figure 3: Self-similar reflection of a straight vertical shock in a convex corner. 
Different "?" patterns occur depending on corner angle and other parameters. 

model equation. Moreover, self-similar flow patterns occur naturally in various 
reflection experiments. 

In self-similar flow, density and velocity are functions of ^ = x/t and rj = y/t 
rather than x,y. To produce a reflection, we consider the horizontal upstream 
wall A and the downstream wall B (see Figure |3|), meeting in the origin and 
enclosing an angle 180° — 9. For t < a vertical incident shock approaches the 
corner from the left, reaching it at i = 0; for t > it continues along B, while 
a complex pattern is reflected back from the corner. For regular reflection, the 
incident and reflected shock meet in a point ^. An observer travelling in the 
reflection point will observe a flow expanding at a constant rate, approaching a 
local RR as in Figure [T] left as t 1 +oo. 

To understand self-similarity intuitively, focus on the corner between the two 
walls in Figure [3] right, t "f oo corresponds to zooming into the corner whereas 
t I corresponds to zooming infinitely far away from the corner. 

The three transition criteria discussed for steady reflection specify angles 9^ 
(detachment), 9s (sonic) and 6'jv (von Neumann), depending on 7 and Lq, so 
that RR occurs for larger 9 whereas IRR occurs for smaller 9. (Here, Lq is 
the Q region Mach number as seen by an observer traveling in the intersection 
point of incident shock and B (= reflection point, in the RR case); of course 
an observer stationary in the corner will perceive a different velocity in the Q 
region.) Note that 6^ < 9s, 9^ always. Figure [H] compares the criteria in the 
case of monatomic gas (7 ~ 5/3). 

It has also been proposed that the correct criterion may not be the same in 
steady and self-similar flow (see below), or tha t there m ay be bistable cases 
where RR and IRR can both occur (see |HOS79lllVF+0l] ). 
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Figure 4: Left: regular reflection. The dotted arc separates a region of constant 
velocity (above) from a nontrivial region. Self-similar potential flow changes 
type from hyperbolic (above) to parabolic to elliptic across the arc. Right: 
single Mach reflection. 



gamma=5/3 reflection (full Euler) 



gamma=5/3 reflection detail (full Euler) 
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Figure 5: Left: transition angles predicted by each criterion (sonic and detach- 
ment almost coincide); right: detail. 
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Nevertheless, it seems that there is an overall preference for the sonic criterion 
in the scientific community, at least for self-similar reflection. 

Numerical and physical experiments are hampered by various difficulties and 
have not been able to select the correct criterion. For example numerical dissi- 
pation or physical viscosity smear the shocks and cause boundary layers that in- 
teract with the refiection pattern and can cause "spurious Mach stems" jWC84| . 
Moreover, 9d and 9s are only fractions of a degree apart (see Figure [S] right), a 
resolution that even sophisticated experiments (e.g. [LD89| ) have been unable to 
reach. To quote |BD92| : "For this reason it is almost impossible to distingiush 
experimentally between the sonic and detachment criteria." 

Constructing exact solutions of most genuinely multi-dimensional fiow problems 
is infeasible or restricted to severely simplified equations. Moreover it would be 
prohibitively expensive if it could only confirm results that have already been 
obtained many orders of magnitude faster by numerical or physical experiments, 
unless the certainty of mathematical proof is needed. Regular reflection appears 
to be the first instance where rigorous analysis might make a genuine contribu- 
tion by answering a problem that could not be resolved unambiguously by other 
techniques. 



1.3 Results 

In this article, using techniques developed in |EL07| . regular reflection is con- 
structed as a self-similar solution of compressible potential flow, with polytropic 
(7- law) gas. While classical rcgular/Mach reflection studies vertical incident 
shocks, wc consider the non-vertical cases too (these may not arise from any 
t < flow), including cases where ^ > f • 

Most importantly, for some values of 7 and upstream Mach number Mj, in 
particular 7 = 5/3 and Mj = 1, every 6 near Og can be covered. This shows 
rigorously that criteria stronger than the sonic criterion are false, at least for 
potential flow with this choice of parameters. 

As discussed above, there is some tendency to believe that regular reflection 
does not persist beyond the sonic criterion; ongoing work aims to show this 
rigorously, at least under mild assumptions. This would rule out the weaker 
criteria as well, in particular the detachment criterion, hence prove that sonic 
is correct. The problem of weak vs. strong reflection (see above) would vanish 
as well. 

However, for now the success is qualified: potential flow lacks contact discontinu- 
ities, so after the transition to (say) SMR the flow pattern must be qualitatively 
different from the full Eulcr flow. It is still possible that the two models may 
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have different transition criteria (however, the author beheves that this is not 
the case). 

Although some genuinely multi-dimensional exact solutions have been con- 
structed for steady Euler flow, self-similar Euler flow is an open and inherently 
rather difficult problem. But again, it seems unlikely that numerical or experi- 
mental techniques will yield a clear — let alone universally accepted — answer 
soon, so rigorous analysis would be very valuable. 

Here is the precise result: 

Theorem 1. Consider potential flow, as discussed in Section \1.5\ Consider a 
wall A ~ (— oo, 0) X {0} (see Figure\^, a second wall ray B at a clockwise angle 
180° — 6 from A, and an incident shock Q, at a clockwise angle 180° — /3q from 
A, meeting B in the reflection point ^r. Assume that there is a corresponding 
reflected shock R in ^/j, emanating down and left (or vertically down). Define 

Vi := {{x,y) e IR^ : y > 0, -cxX a: < ycotiPq)} 

Vq := {{x,y) e : y > 0, ycot(/3Q) <x< ycotO}, 

V := {(x, y) G IR^ : y > 0, -oo <x <ycotd} 

(see Figure^. 



1. Assume the sonic criterion holds: > 1 in in the sector below R. 

2. Assume that 

\vi • 7is| < CI (1.3.1) 

3. Envelope condition: of the tw^ points on the R shock with Ld = 1, let 

be the one closer to ^r. Consider shocks with upstream data vi,pi that 
go from counterclockwise and satisfy Ld < 1 in every point. Assume 
that all such shocks reach A before meeting B or the circle with center vj 
and radius cj. 

Then there exists a wea solution = c^{t,x,y) G C"'H[0,oo) x V) of 

unsteady potential flow for t > 0, x € V , (1.3.2) 

• n = on dV, (1.3.3) 

p = Pu V(i) = vi fort^O,xeVj, (1.3.4) 

p = PQ, S7cb = VQ fort^O, xeVQ. (1.3.5) 



^see Section ll.6l 
^see Remark ll.3.11 
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Figure 6: Left: a local RR pattern; right: the curved portion S of the reflected 
shock has Ld < 1, hence must be left of the envelope E, which bounds it away 
from the dotted circle and from B. 
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Figure 7: The initial data is constant in each of two sectors that are separated 
by the incident shock Q 
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Of course existence by itself merely validates that potential flow has interesting 
solutions. In addition, detailed results about the structure of the weak solution 
can be obtained (see Rcmark l2.9.ip : most importantly, the flow patterns are of 
RR type. 

Remark 1.3.1. By weak solution we mean that 



for all test functions ■& G C^i^). 

(For (j) g C"'^(ri), the velocity V0 is a.e. well-defined on {0} x V, but 0t and 
hence p may not be well-defined.) 

Remark 1.3.2. Condition (|1.3.ip and the envelope condition arc merely techni- 
cal. The envelope condition is needed in some cases to prove the shock docs 
not vanish (which is never observed in numerics); none of the other estimates 
requires it. Both conditions can probably be removed by future research. 



1.4 Related work on constructing exact solutions 

In recent years multi-dimensional compressible inviscid flow has received re- 
newed attention, after several recent breakthroughs brought the theory of one- 
dimensional compressible flow to a satisfactory state jGli651[BB051ILY99llBCP00| . 

|EL07| (see also jEL05bl[EL06| ) studies supersonic flow onto a solid wedge. For 
sufficiently sharp wedges, the steady solution consists of a straight shock on each 
side of the wedge, emanating downstream and separating two constant-state 
regions. In inviscid models this shock wave must keep the downstream velocity 
tangential to the wedge surface (slip condition). As for regular reflection, there 
are two different shocks for each (small) wedge angle, a weak and a strong shock. 
The weak shock is more commonly observed, but no mathematical argument 
was known to favor it prior to [EL 07] . In that article, an exact solution was 
constructed for a wedge at rest in stagnant air, accelerated instantaneously to 
(sufficiently high) supersonic speed at time 0. The resulting flow pattern is 
self-similar and has a weak shock at the wedge tip. 

Many of the techniques in [EL 07] are essential in the present article. 

The most closely related work, and so far the only other paper that proves global 
existence of some nontrivial time-dependent solution of potential flow is [CFj : 
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V0(O, x) = VI 
V0(O,f) = i7q 



for a.e. x G Vj 
for a.e. x E Vq 



(1.3.6) 
(1.3.7) 



and 
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using different tecliniques, they construct exact solutions for regular reflection, 
assuming sufficiently blunt wedges (6* « ^). 

Some prior wor k studies refl ection and other problems for simplified models of 
gas dynamics. |CanicKK02j consider regular reflection for the unsteady tran- 
sonic small disturbance equation as model. [Zhe06| studies the same problem 
for the pressure-gradient system. The monographs [ZheOli ILZY98| compute 
various self-similar flows numerically and present some analysis and simplified 
models. 



1.5 Potential flow 



Here we briefly present derivation and elementary results for potential flow. 
More information can be found in [EL 07] . 

Consider the iscntropic Eulcr equations of compressible gas dynamics in d space 
dimensions: 

pt + \/-ipv)^0 (1.5.1) 

d 

{pv)t + Y^ipv'v),. + \/{p{p)) = 0, (1.5.2) 

i=l 

Hereafter, V denotes the gradient with respect either to the space coordinates 
X = {x^,x^, • • • , x'^) or the similarity coordinates t~^x. v = {v^ ,v'^ , ■ ■ ■ , v"^) is 
the velocity of the gas, p the density. p{p) pressure. In this article we consider 
only polytropic pressure laws (7-laws) with 7 > 1: 



7 \Po, 

(here cq is the sound speed at density po). 

For smooth solutions, substituting (|1.5.1|) into (|1.5.2p yields the simpler form 

vt+v-\/'^v + \/{n{p))^0. (1.5.4) 

Here tt is defined as 

-, 7 > 1 



nip) = Co . 



7-1 

log(p/po), 7=1. 



This TT is C°° in /9 e (0, 00) and 7 G [1, 00) and has the property 
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If we assume irrotationality 



= vi 



(where i, j = 1, . . . , d), then the Euler equations are reduced to potential flow: 

V = Vgcj) 

for some scalar potential function cj). For smooth flows, substituting this into 
(|1.5.4p yields, for z = 1, . . . , d, 

= 4t + V4 • V(/) + ^(p), = (0t + + ^(p))^. 

Thus, for some constant A, 

p = 7r-i(A-0,-^). (1.5.5) 

Substituting this into (|1.5.ip yields a single second-order quasilinear hyperbolic 
equation, the potential flow equation, for a scalar held <f>: 

{pi^t, |V0|))^ + V • (p(0t, \V^\)V4>) = 0. (1.5.6) 

Henceforth we omit the arguments of p. Moreover we eliminate A with the 
substitution 

^ ^ 0, (/)(i, x) <- (i){t, x) - tA 
(so that (j)t ^ (pt — A) . Hence we use 

p = ^-i(-0t-i|V0n (1.5.7) 

from now on. 
Using = Pp and 

(.-!)' = (.,)-! = (^)-i = 4 (1.5.8) 
the equation can also be written in nondivergence form: 

d 

(Ptt + 2\/(j>t • V0 + ^ MJ<I>^J - c^^4> = (1.5.9) 

(|1.5.9[) is hyperbolic (as long as c > 0). For polytropic pressure law the local 
sound speed c is given by 



c' = c^ + (7-l)(-</'t--|V0n. (1.5.10) 



*We consider simply connected domains; otherwise might be multivalued. 



February 5, 2008 



12 



1.5 Potential Row 



1 INTRODUCTION 



Our initial data is self-similar: it is constant along rays emanating from x = 
(0,0). Om' domain V is self-similar too: it is a union of rays emanating from 
{t,x,y) ~ (0,0,0). In any such situation it is expected — and confirmed by 
numerical results — that the solution is self-similar as well, i.e. that p, v are 
constant along rays x = t^ emanating from the origin. Self-similarity corre- 
sponds to the ansatz 

(j){t,x) ■.= t^jj{i), ^:^t-^x. (1.5.11) 
Clearly, e C°'^{n) if and only if V e C"'^iCW). This choice yields 
v{t,x) = V(t){t,x) ^ Vipit-^x), 

Pit, x) = 'K-\-4>t - i|v0n = ^-i(-^ + - i|vv^n. 

The expression for p can be made more pleasant (and independent of ^) by using 
this yields 

P = ^-'(-X-^|Vxn. (1.5.12) 
Vx = VtA — C is called pseudo-velocity. 
(|1.5.6p then reduces to 

V-(pVx) + 2p = (1.5.13) 

(or -\-dp, in d dimensions) which holds in a distributional sense. For smooth 
solutions we obtain the non-divergence form 

[c^I - VxVx^) : V^x = (c' - X|)X« " 2x«Xr,X?, + (c' - X^Xm = |Vxl' - 20^ 

(1.5.14) 

Another convenient form is 

[c^I - VxVx^) : V^^ = (c2 - xD^a - '^XiXni'in + (c' ^ x'^w = 0. 

(1.5.15) 

Here, (|1.5.10p for polytropic pressure law yields 

c^=cl + {^-l){-X-\\Vx?) (1.5.16) 
Remark 1.5.1. ()1.5.13p inherits a number of symmetries from (jl.S.ip . (|1.5.2p : 

1. It is invariant under rotation. 
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2. It is invariant under reflection. 

3. It is invariant under translation in ^, which is not as trivial as translation in 
x: it corresponds to the Galilean transformation w <— u + vq, x <— x — vqI 
(with constant vq £ IR'') in {t, x) coordinates. This is sometimes called 
change of inertial frame. 

(|1.5.14p is a PDE of mixed type. The type is determined by the (local) pseudo- 
Mach number 

L:=^, (1.5.17) 
c 

with < L < 1 for elliptic (pseudo-subsonic), L = 1 for parabolic (pseudo- 
sonic), L > 1 for hyperbolic (pseudo-supersonic) regions. 

While velocity v is motion relative to space coordinates x, pseudo-velocity 

z:= Vx 

is motion relative to similarity coordinates ^ at time t = 1. 

The simplest class of solutions of (jl.5.14p are the constant- state solutions: ip 
affine in ^, hence v, p and c constant. They are elliptic in a circle centered in 
^ = V with radius c, parabolic on the boundary of that circle and hyperbolic 
outside. 

Convention 1.5.2. If we study a function called (e.g.) x, then ip, p, L etc. 
will refer to the quantities computed from it as "01 Pi are computed from x 
(e.g. "0 = X + ^ICP)- We will tacitly use this notation from now on. 



1.6 Potential flow shocks 

Consider a ball U and a simple smooth curve S so that U = J7" U U J/'^ where 
U^,U'^ are open, connected, and S,U'^,U'^ disjoint. Consider x ■ U \R so 
that X = X"-'^ in where x"''^ e C^{IP^). 

X is a weak solution of p.5.13p if and only if it is a strong solution in each point 
of U- and U+ and if it satisfies the following conditions in each point of S: 

X"=x', (1-6.1) 
n • (p^Vx" - p'^Vx'') = (1.6.2) 

Here n is a normal to S. 
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(|1.6.ip and (|1.6.2p are the Rankine-Hugoniot conditions for self-similar potential 
flow shocks. They do not depend on ^ or on the shock speed explicitly; these 
quantities are hidden by the use of x rather than ip. The Rankine-Hugoniot 
conditions are derived in the same way as those for the full Euler equations (see 
|Eva98[ Section 3.4.1]). 

Note that (|1.6.ip is equivalent to 

tA"* = i'''. (1.6.3) 
Taking the tangential derivative of (|1.6.1|) resp. (I1.6.3P yields 



dx" dx'^ 
dip"" _ dip"^ 



(1.6.4) 
(1.6.5) 



The shock relations imply that the tangential velocity is continuous across 
shocks. 

Define (zjj, z^) := := Vx" and (u^, vD := Vu := V^"^. Abbreviate := z„ -t, 
:= Zu ■ n, and same for v instead of z. Same definitions for d instead of u. 
We can restate the shock relations as 

p„z: = pdz2, (1.6.6) 
4 = 4- (1-6.7) 

Using the last relation, we often write 2* without distinction. 

The shock speed is tr = ^ • n, where ^ is any point on the shock. A shock is 
steady in a point if its tangent passes through the origin. We can restate ()1.6.6p 
as 

- Pdv2 = a{pu - pd) 

which is a more familiar form. 

We focus on pu, pii > from now on, which will be the case in all circumstances. 
If pu ^ Pd in a. point, we say the shock vanishes; in this case z^ = z^ in that 
point, by (|1.6.7|) . In all other cases z^,z^ must have equal sign by p.6.7p : we 
fix n so that z^, z"^ > 0. This means the normal points downstream. The shock 
is admissible if and only if p„ < pd which is equivalent to z^ > z^. 

A shock is called pseudo-normal in a point ^ if z* = there. For — 0, this 
means that the shock is normal (u* = 0), but for ^ 7^ normal and pseudo- 
normal are not always equivalent. 
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It is good to keep in mind that for a straight shock, pd and Vd are constant if 
Pu and Vu are. Obviously Zd may vary in this case. 

We will need two detailed results. 

Proposition 1.6.1. Consider a fixed point on a shock with upstream density 
Pu and pseudo-velocity Zu held fixed while we vary the normal. Define (3 := 
Z(z„,n). pd is strictly decreasing in \f3\, whereas Ld, \zd\ are strictly increasing. 
Cd is strictly decreasing for "f > 1, constant otherwise. Moreover 

{di,Vd) ■ ri - {d^zd) . = - l), (1.6.8) 

(dpvd) ■ t = (dpzd) ■ r= z^^ - zl. (1.6.9) 

If Zu = {z^^Q) with z^ > 0, then z^ is increasing in 

Proof This is |EL07[ Proposition 2.5.1]. □ 

Proposition 1.6.2. Consider a straight shock with ~ 0, v-^ < and down- 
stream normal n = (sin/3, — cos/3) through ^ = {0,ri).For every f3 G (— ^,-|) 
there is a unique rj = r]^ €z IR so that = 0. 77Q and the corresponding down- 
stream data are analytic functions of p. rjQ is strictly increasing in \(3\. 

For the shock passing through (0, 77q), let ^2 and be the two points with Ld = 
•\/l — £■ These points are analytic functions of (3. L", pd and z" are increasing 
function^ of (3; and LJJ are decreasing functions of p. For (3 £ [0, -1), vfl 
is a strictly decreasing function of (3 with range (?7* , 77o], where tJq is the r]^ for 
P = 0, and Tj*^ is some negative constant. 

Proof. This is |EL07|. Proposition 2.6.2]. □ 
1.7 Envelope 

Many techniques in this paper are similar to the construction in [EL07| ; Section 
4.2 in loc.cit. is a good overview. However, in [ELOTi Proposition 4.11.1], a 
lower bound for the shock strength is obtained by a delicate argument using the 
density. Although this argument would reproduce the results of |CFj (namely 
RR existence for 6 « ^), it cannot prove the main new contribution of this 
paper: existence (at least in some cases like 7 = 5/3, M/ = 1) of RR for 6 k, 0^ 
(with > 6s), where 9s is the smallest 9 allowed by the sonic criterion (see 
Section [L^ . 

^AU of these are independent of the location along the (straight) shock. 
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Figure 8: Left: through each ^ farther than c„ from t/„ there are exactly two 
straight shocks (solid, dashed) with Ld = I, mirror images of each other. The 
shocks with Ld < 1 are between them (indicated by arc left of i). The sohd 
lines define the direction field whose integral curves are "envelopes". Right: no 
shock with Ld < 1 can approach w„ faster (in counterclockwise direction) than 
the counterclockwise envelope. 

For this goal, a new idea is needed: as we will show, the curved portion S of 
the reflected shock in Figure S] left has an elliptic region of potential flow on 
its right (downstream) side, hence downstream pscudo-Mach number Ld < 1 
everywhere. Such a shock cannot vanish until it reaches the circle of radius c/ 
around vj; moreover < 1 is a constraint on the possible shock tangents, so 
that the shock cannot reach the circle quickly. It is bounded away from the 
circle by the envelope: 

Definition 1.7.1. Given constant upstream velocity Vu and sound speed €„. 
Consider a shock through a point ^ with = — ^| > c^. As shown in 
Proposition II. 6. n Ld is strictly increasing in \(3\ where (3 = Z(£u,n) G (— 7r,7r] 
is the counterclockwise angle from z„ to n. 

There are exactly two shock normals so that Ld = I- They arc mirror- images of 
each other under reflection across the line with tangent z„ through ^ (see Figure 
[8] left). Consider the one with /? > 0; its tangent spans the solid line on Figure 
[5] left. The tangents for different ^ form a direction field. The counterclockwise 
envelope is defined to be a maximal integral curve of that direction field (see 
Figure [8] right) . 

We can parametrize the envelope (like other smooth shocks) in polar coordinates 
(r, (j)) centered in ?;„ , by a function (p i-^ r* ((/>) (because the shock relations do 
not admit shocks with a tangent passing through z4). The counterclockwise 
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envelope satisfies an ODE of the form 

= (1.7.1) 

for some analytic /. 

Wc will not need the fact, but explicit formulas for / can be derived. For 
example for 7 > 1, 



/(••) = 



1 7+1 ( 7+1 

^ 7-l+2(r/c„)-^ 1^2+(7-l)(r/c„)^ 
7 + 1 



(1.7.2) 



7-l+2(r/c„)' 



Moreover it can be shown that the envelope always reaches the circle, meeting 
it in a point where the envelope is C^, but not more regular, and tangent to the 
circle; it cannot be continued beyond that point. 

Proposition 1.7.2. Let some smooth shock be parametrized as (f> s- r(t/)); let 
the envelope be parametrized by cj)^ r*((j)). Assume that < 1 in every point 
of the shock. If r{(j)o) > r*{4>o) for some 4>q, then r{(j>) > r*{(j)) for (j> > (f>Q. If 
instead Ld > 1 in every point of the shock, then r{4>) < r*{4>) for 4> > 4>o. 



Proof. Our discussion above can be restated as follows: Ld < 1 for the shock 
means — /3* < (3 < (3* where (3* is the (3 for the envelope. Hence 

\%\<nrm. 

In particular 

dr 

^ > -firm. 

Since / is smooth, in particular Lipschitz, the invariant region theorem shows 
that the shock cannot meet the envelope for > c/jq. □ 



In Proposition 12.6.21 we will exploit this fact to bound the curved portion of 
the reflected shock away from the downstream wall and to ensure its uniform 
strength. 



1.8 Sonic criterion 



We focus on the classical case of vertical incident shocks. In some cases. Theorem 
[1] allows us to construct a regular reflection pattern like Figure |4] left for every 
myrefsection:refl). As 6 [Og, the dotted6' > 6s near 9s, where Os is the smallest 
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Figure 9: For (3q = (vertical incident shocks) and the set of M/,7 enclosed 
below the dashed and above the solid line, solutions can be constructed for all 



allowed by the sonic criterion (see Section parabolic arc in Figure |4] left 
approaches the reflection point. 

To check whether the envelope condition is satisfied for a particular choice of 
9 and incident shock, it suffices to find the reflected shock and on it (see 
Theorem [1]) and to integrate the ODE p.7.ip defining the envelope. Although 
the ODE is trivially separable, the resulting integral and nonlinear algebraic 
equation do not have an explicit solution except for special values of 7 (see 
(|1.7.2[) '). Numerical integration is needed to check whether the envelope meets 
B or the circle with center vj and radius c/ before it meets A. 

In Figure [H we consider arbitrary 7 e [l,oo) and A// <E (0, 00) while fixing 
9 = 9s- Values of 7 and Mj above the dashed curve do not admit a vertical 
incident shock with zero velocity in the Q region (a similar phenomenon occurs in 
the full Euler equations). Values below both solid and dashed curve violate the 
envelope condtion. Values between solid and dashed curve do have an incident 
shock as well as a reflected shock that satisfies the envelope condition. 

The smallest possible 7 in that feasible region is 7 = 1.626354... with Mj = 1. 
In particular the monatomic gas case 7 = 5/3 is covered, whereas 7 = 7/5 or 
7 = 4/3 arc not covered. (However, the latter values are also possible if we allow 
non-vertical incident shocks.) For 7 = 5/3, Mj = 1 we have 0s = 55.4583...°; 
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for 9 = 9s the envefope meets A in the point (—0.000012..., 0), just enough to 
avoid B and the circle. 

While the proof of Theorem [T] itself is rigorous, checking the envelope condition 
is done numerically here, i.e. not a mathematical proof in the strict sense. How- 
ever, the shock relations form a small system of nonlinear algebraic equations 
and the envelope is defined by (|1.7.ip . a scalar nonlinear ODE which is benign 
except for a mild singularity as r | 1. The numerical methods for these types 
of equations are well-understood and a complete convergence theory and error 
analysis is available — which is not at all the case for the full Euler or poten- 
tial flow PDE. Another option is to study rigorous proofs in various asymptotic 
limits such as Mj | 0, 7 t oo. Moreover the envelope condition is most likely 
unnecessary since regular reflection up to = 9^ is observed in numerics for 
many other values of 7 and Mi as well. Since we expect that the condition 
will be eliminated by further research, it makes little sense to strive for absolute 
rigour at this point. 



2 Construction of the flow 

The elliptic region is constructed as follows: we define a function set T by im- 
posing many constraints on a weighted Holder space C^'" (weighted to account 

for loss of regularity in the corners). An iteration /C : — *■ C^'" is constructed 
so that its fixed points solve the PDE and boundary conditions for the elliptic 
region (see Remark |2. 2. 4p . !F and K. depend on several parameters like 7, col- 
lected in a parameter vector A. To show that /C has a fixed point for all A, we 
use Leray-Schauder degree theory. 

Most of the effort is spent on showing that K. does not have fixed points on 
dT, which implies that K, has the same Leray-Schauder degree for all A. As 
dJ-^ is defined by constraints in the form of inequalities with continuous sides, 
this is achieved by showing that a fixed point satisfies the strict version of each 
inequality (< instead of <). 

A major technical difficulty are the parabolic arcs (dotted arc in Figure [4] left) 
where self-similar potential flow ()1.5.15p degenerates from elliptic to parabolic. 
This problem has been solved in |EL07| (and, by different techniques, in |CF| ). 
by modifying the arc to be slightly elliptic, with boundary condition = 1 — e, 
and obtaining estimates uniform in e. 

For a particular choice of A the problem is much simpler (see Figure [T5|) . In that 
case an explicit solution can be given and shown to be unique and have nonzero 
Leray-Schauder index. This implies that K. has nonzero degree, hence at least 
one flxed point, for every A. The flxcd point is extended to a solution on the 
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VI 




Figure 10: Perturbation from the trivial case of R parallel to the wall. 



entire domain by adding the hyperbolic regions and interface shocks. Using the 
e-uniform estimates as well as compactness, we can pass to the limit e J, to 
obtain a solution of our problem. 



2.1 Parameter set and definitions 

Instead of working in the setting of Theorem [U it will be convenient to choose 
parameters in a different way. 

Choose pi,ci > 0. Note that we may fix pq and Cq in the pressure law (|1.5.3p 
separately; however, given these constants (and 7), every other c is a function 
of p only (and vice versa). 

Let e > be sufficiently small for the following. Consider a vertical downward 
velocity vj onto a solid wall B (see Figure [T0|). According to Proposition 11.6.2] 
there is exactly one straight shock with upstream velocity Vu ~ vi and sound 
speed Cu = cj so that Vd = 0; that shock is horizontal. Let 77^ > be its 
vertical coordinate. Of the two points on that shock with Ld = \/l~~e, let 
^ = (Ccj^c) t^^*^ right one. By the same proposition, the shock belongs to 
a smooth one-parameter family of shocks, each called R shock, parametrized by 
7?J e (0, 1]^], so that =0 and so that = (Cc' Vc) ^^e right Ld = e 
point. Define Mj E [—1,0) to be w|/c/ in these coordinates. Note that (|1.3.ip 
rules out Mj < — 1. Let vn ~ Vd be the downstream velocity of the R shock. 

It is not clear whether there is an incident shock Q matching each refiected 
shock R. In fact for tjq = 77^, the R shock does not even meet B_, so clearly 
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Figure 12: To avoid degeneracy, we impose a "slightly elliptic" boundary condi- 
tion, = 1 — e for e > 0, on P'*^' . The shock S is free, along with the endpoints 
and which may slide freely on A resp. P(^). But for fixed points can 
be shown to be close to hence to ^*(°\ 
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there is no RR. However, for the construction of the elliptic region, a Q shock 
or reflection pattern are not needed. 

To complete the situation of Theorem [1] a wall A is needed. To satisfy the slip 
boundary condition (u/ — ^) • n = on A, necessarily the extension of A to a 
line has to pass through vj. We fix A by choosing ^ab on B. 

Let E be the counterclockwise envelope starting in ^c- If ^ meets 5 before it 
meets the circle with center vj and radius c/ (Figure [TT] left), let ^eb be that 
point. Otherwise (Figure [Tl] right) take the line through vi and the meeting 
point of E and circle, and let £^eb be its intersection with _B. We allow 

UBei^EB^V^B] (2.1.1) 

(and r]AB = obviously). This constraint ensures that (1) the envelope meets 
A first, while (2) R and A form a sharp or right angle. 

Given ^ab we let B be the part of B right of £,ab ■ A is the half- line upwards 
starting in ^ab whose extension passes through vj. Let Ha be the unit normal 
of A pointing left, Hb the unit normal of B pointing down. Let be the 
downstream (hence downwards) unit normal of the R shock. For each n?, t-? is 
always the corresponding unit tangent in counterclockwise direction. 

Remark 2.1.1. Every local RR pattern that satisfies the conditions of Theorem 
[T]is covered by the parameter ranges defined above. 

Pi and vj define a potential ^j/ for the / region: 

Similar potentials -0^ and tp'^ (if an incident shock Q exists) are defined by 
Pfl, Wfl and pq,vq. 

Now we use Remark 11.5.11 invariance under translation. Translation in self- 
similar coordinates corresponds to a change of inertial frame, i.e. to adding a 
constant velocity to all u, ^. Moreover we may rotate by Galilean invariance. 
This changes Figure [TT] to Figure [12] which has the coordinates in which we 
originally posed the self-similar reflection problem. 

Let P^"^ be the circle arc centered in vn with radius cji ■ \/l — e (see Figure 
1121 where the coordinates have been changed), passing from on B counter- 
clockwise to ^"^"^ on R, excluding the endpoints. (Wc omit the superscript e if 
it is clear from the context.) will be called the expected corner location. Let 
B'^'^^ be the part of B from (^ab to (excluding the endpoints). 

Take n/j, uq to be the downstream unit normals of the shocks R, Q (nj^ points 
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towards B). Let nA,nB be outer unit normals of A,B, i.e. pointing away from 
the gas-filled sector V enclosed by B,A. 

We choose an extended arc P that overshoots by an angle Sp > 0, which 
we choose continuous in 7,Ca_b,??c- The particular 5p is not important, but it 
may not depend on e, and P may not have a horizontal tangent in Figure [Ml 
coordinates. 

P*, P, and later P, are called quasi-parabolic arc (or parabolic arcs, by abuse 
of terminology, or short arcs). 



Parameter set The Definitions I2.1.2[ 12.2.21 and 12.2.31 use many constants 
and other objects that will be fixed later on. In all of these cases, an upper 
(or lower) bound for each constant is found. Whenever we say "for sufficiently 
small constants" (etc.), we mean that bounds for them are adjusted. To avoid 
circularity, it is necessary to specify which bounds may depend on the values 
of which other bounds. In the following list, bounds on a constant may only 
depend on bounds on preceding constants. 

Sp, Cl, CrpSsBi Scc, Spc,5pn,5d, Sp, SLb, 

Cpt, CytR, CynA, Gsn, SvtA, SvnB,So, Cd, £, Cq, f/, a, /?. (2-1-2) 

The constants Cc,ri,a, (3 may depend on e itself, not just on an upper bound. 
ri may also depend on "0. The reader may convince himself that the remainder 
of the paper does respect this order. 

The parameters 7, and ^ab used in Leray-Schauder degree arguments will 
be restricted to compact sets below so that any constant that can be chosen 
continuous in them might as well be taken independent of them. Dependence 
on other parameters like pi will not be pointed out explicitly. 

Constants S-? as well as a, (3,ri,e are meant to be small and positive, constants 
C? are meant to be large and finite. 

Definition 2.1.2. For the purposes of degree theory we define a restricted 
parameter set 

A:= {a= (7,77c,Cas) : 7 S [1, 7], S [VcXcIUb £ Hab^Iab]} 

where it is important that ^ab and vf^ are the values in the coordinates of 
Figure [TO] and Figure [TTl clearly their values are entirely different in any other 
coordinate system we use. 7 € [l,oo) is an arbitrary constant. Moreover, 

— * I 0, 7 = 1: /r, 1 o^ 
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sector of 




■CI 

Figure 13: The shaded sector consists of all A rays that (1) form a sharp angle 
with R, while (2) meeting E before E meets _B or the circle. 



and 

where Cj, Cjj (to be determined in Proposition 1 2. 5.4]) do not depend on e or A. 

is a constant satisfying < ry^ < Ty^. Finally, € (C-EB, Cab] depend 
on 7 and 77 J. 



Proposition 2.1.3. A contains {^,7]^, S^ab) = (Ij '/ci '^^li) and is path- connected, 
for e sufficiently small (depending on Crj, C^) and ^ sufficiently close to £,eb- 



Proof. 1. Wc note that the interval {£,eb,Vij] has boundaries that are con- 
tinuous functions of A. 

2. The interval is always nonempty: consider the coordinate system and 
setting of Figure [Til extended in Figure [131 Consider a line T through 
that touches the (upper half of the) circle with center z// and radius 
c/ in a point ^y. T can be considered a zero-strength shock (velocity iji, 
density pi on both sides), with = 1 in and L^> 1 elsewhere. Hence 
Proposition 11.7.^ applies : let (/) i~> r(0) parametrize the line segment from 

^'^^ to ^T] let (j) I— *■ rE{4i) parametrize E. Then rE{4>) > r{4>) > c/ on the 
interior of the corresponding (j) interval, so the envelope E cannot touch 
the circle right of <^t- Moreover, since we have assumed that Mf < 1 
(restriction (|1.3.ip ). that means the circle cither meets or intersects B. If 
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E meets B_ before it meets the circle, then necessarily it meets the part of 
B left of the circle first. 

On the other hand, the extremal choice ^ab = vf, for A corresponds to (a 
segment of) the line through vi and vr (right side of the shaded sector in 
Figure I13[) , which is perpendicular to R. Its intersection with the circle is 
necessarily right of ^t- Thus: if E meets the circle before it meets B_, then 
S,EB < necessarily. If E meets S before the circle, then it must meet 
it left of the origin, so £,eb < < v^. Either way the interval {£,eb, vf/] is 
nonempty. 

Threfore the interval vf/^d^-e^^^] is also nonempty, if e is sufficiently 

small (depending on Cj) and sufficiently close to ^eb- 

3. Finally, we show that the special A = (1, rj^, vf^) can be connected by paths 
in A to all other A: it connects to any (1,77^,^^5) with rj^ £ iv'^iVc) 
and Ub G Hab^Vr]- These include (l,r;[l, ~ C,, • - Q~- e'^') 

which connects to any (7,??^ — C^, • e^/^,u^ — C^e^/^) with 7 > 1. This 
point, in turn, connects to any , i]^ , £,ab) with 77^ G [tIqiVc] and £,ab € 
[£,^g,^AB]- Hence A is path-connected. 

□ 



2.2 Function set and iteration 



Definition 2.2.1. Let U C IR" open nonempty bounded with dU uniformly 
Lipschitz. Let F C dU. For k E INq, a S [0, 1] and /3 G (—00, k + a] we define 
the weighted Holder space C^'"(f7, F) as the set of m G C'='"(f7 - F) so that 



is finite. 



Definition 2.2.2. For sufficiently small dp > 0, there is a function b e C^{V) 
with b, |V6| < 1 so that 6 = on P*"', 6 > elsewhere, 6„ on 1 and B, 
and so that b depends continuously on the parameters A but is independent of 
e. From now on we fix a particular b. 

Proof. The construction is straightforward. 6p is taken so small that P^^^ does 
not meet Au BU {^ab} except in . □ 

Definition 2.2.3. 
Onion coordinates 



February 5, 2008 



26 



2.2 Function set and iteration 



2 CONSTRUCTION OF THE FLOW 




VC 



A 




Ub^ {0,0) 



Figure 14: Transformation to "onion" coordinates {a, rf) 

Rotate Figure [T^ so that B is the positive horizontal axis (see Figure [T3] left), 
then shift horizontally so that u/ is vertical (Rcmark ll.S.ip . Define new coordi- 
nates (cr, 77) e IR^ (see Figure [T3] right) so that 

1. the coordinate change from (^,77) to {a,ri) is C°° with C°° inverse, 

2. B(<^) maps to (0, 1) x {0}, 

3. A maps to {0} x (0,oo), 

4. P^'^^ maps to {1} x (0,?7c) (where fjc is the -q coordinate of the upper 
endpoint of P'*^'), 

5. maps to (1, 0), 

6. S^AB maps to (0,0). 



We require that the change of coordinates and its inverse depend continuously 
(in the C°° topology) on A G A. The construction is straightforward. 

Here and in what follows, we will use the weighted Holder spaces C^'"(f/), as 
in Definition HXlJ The domain U is either [0, 1]^ with F = {(0, 0), (1, 1)}, or 
with F ~ {£,AB,£,c} (to be defined). For the shock parametrization we use 
U = [0, 1] with F = {0, 1}, or (in Figure [Til left coordinates) U = [Ca,Cc] with 
F = {^c}; for f/ = P we use F = {^c}, and for U = A oi U — B wc take 
F = {£,ab}- We omit F as it will be clear from the context. (3 e (1,2) and 
a G (0, /3 — 1] will be determined later. C^'" are Banach spaces so that standard 
functional analysis applies. Moreover, C^'"(r2) is continuously embedded in 
C^{VL), so we have regularity in the corners as well, which is crucial. 
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Free boundary fit 

Let T be the set of functions ip e C^'"([0, 1]^) that satisfy aU of the many 
conditions explained below. Require 

mcl-i[o,m<Cc{^)- (2-2.1) 

The curves of constant a (isolines) in the 77) coordinate plane are nowhere 
horizontal, since the other coordinate is 77. Moreover tp^ = < and V'| = 
vf = 0, so for all a G [0, 1] there is a unique point s(cr)) on the isoline so that 

V''(^,s(a)) = 7A(a,l). (2.2.2) 

We define another coordinate transform by first mapping (cr, G [0, 1] to {a, rj) 
with 77 = s((t)C and then mapping to ^ with the previous coordinate transform. 

Let resp. be the ^ coordinates for the (cr, C,) plane points (0, 1) and (1, 1). 
Let S be the ^ plane curve for (0, 1) x {1} (it is the graph of s, with endpoints 
£,A and ^c)- Define P resp. A resp. to be the image of {1} x (0,1) resp. 
{0} X (0,1) resp. (0,1) x (0,1). 

Require shock-wall separation: 

diS,B)>SsB>0. (2.2.3) 

(|2.2.3p ensures that the map from (a, () to ^ is a well-defined change of coordi- 
nates, uniformly nondegenerate (depending on SsB and Cc), withC^'"([0,l]2) 

resp. C^'"(0) regularity. It is clear now that dil is the union of the disjoint sets 

S, P, A, B, and {Cc, 63, 61, Cab}. 

Require: corner close to target: 

\ric-vh\<^"^ (2.2.4) 
We require e to be so small that £ P- 

For later use we define 77^ := 77^ ± e^/^ and let be so that G Pc- 
Corner cone: 

sup_Z(e - Cc, f - Cc) < ^ - 5cc- (2.2.5) 

(Z(£, 77) is the counterclockwise angle from x to 77.) 
Iteration 
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Here we change to the coordinates of Figure [T^] for the remainder of the defini- 
tion. 

Shock strength/density: require that 

-X-^|Vxl'>0, (2.2.6) 

so that p is well-defined (see (|1.5.12p ). and require 

imn p > pi + Sp. (2.2.7) 
n 



Pseudo-Mach number bound: require 

<l-6Lb-b in H, (2.2.8) 

(Note that L is well-defined because by (|2.2.7|) p>0, sooO.) 6 = on P^^ 
which has distance > | (for sufficiently small e) from il, so p.2.8p implies 

i'<l-^|V6|L-<5Lh-e<l-i(5L6-e in H, (2.2.9) 



Require: there it[f| a function e C|'"(ri) with the following properties: 

1. ip close to V': 

\H'-M\cl^-aos]^)<riW (2.2.10) 

where r/ G C{J-; (0, oo)) is a continuous function to be determined later. 

2. Right away we require rj to be so small that 

-X-^|VxP>0, (2.2.11) 
so that in particular p is well-defined and positive. Moreover, require 

VtAt^w/, (2.2.12) 

3. We require 77 to be so small that (using (|2.2.9p ) 

(co + (1 - 7)(X + l\^xf))l - Vx' > 0, (2.2.13) 
i.e. is a (symmetric) positive definite matrix. 



^1/1 is the product of an iteration step with input 1/). We will ensure in Proposition I2.2.'7l 
that 1/1 is unique and continuously dependent on ip. 
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4. Let C = C{ijj,t[^) be defined in <^ coordinates as 

((eg + (1 - 7)(X + lm\'))l - Vx') : (2.2.14) 

2 2+(l-e)(7-l) ' 

|u/ - Vip\ 

Vi^-nA,V^-riB). (2.2.17) 



(2.2.15) 



where the codomain is 

Y C^,'^,m X 4_",(:S) X C^-- (P) X C^- (3) X C^^'- (:B). 

p.2.16p is well-defined by (|2.2.1ip and (|2.2.12p . The other components 
have no singularities. 

Note: VtP e so |VxP G C^':!*!, so 

{{c^ + (1 - 7)(X + lm\'))l - ^f) e CL"i C^'^-' 

(a < /3 - 1 as required above), and V^ip € so p.2.14p is e C^-"^^. In 

the same way we check that p.2.15p . p.2.16p and (|2.2.17p are C^f^ 

For 'ip we use the C^'"(ri) topology. We pull back -0 and the value of C to 
(cr, C) coordinates, via the coordinate transform defined by V' (see above), 
so that we have a fixed domain [0, 1]^ for all Banach spaces. Then £ is a 
nonlinear smooth map in the corresponding topologies. 

Most importantly: require 

£(V,V^) = 0. (2.2.18) 

Other bounds 

Require 



lco,i(Ti) < Cl (2.2.19) 



where Cl may not depend on e. 
Xt and Xn on parabolic arc: 



maxc-i|^| <Cpfei/^ (2.2.20) 



p 



dt 
-idx 



maxc~^-^ < -5pn- (2.2.21) 
P on 
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We emphasize that Spt, Spn may depend only on A, but not on e (or ip). 
Velocity components: 

v-nA< CynA ■ e^/^ in H, (2.2.22) 



v-tR< vr ■ tR + C„tR ■ e^/^, in fl, 



(2.2.23) 



vnB 



in 17 



and 



V-tA<Vi - tA - 8ytA hi VL. 



(2.2.24) 



(2.2.25) 



Shock normah Let N C (unit circle) be the set of n counterclockwise from 
Hr to tA- Then the shock normal satisfies 



sup d{n, N) < Csn ■ e 
s 



1/2 



(2.2.26) 



Set El := A, E2 S, := P and S4 
(|2.2.17p . (12.2.16P of/: as 



S. Write the components (|2.2.15p . 



where the ^ dependence includes the dependence on xiO ^nd Vx(C)- 

(7^ has some singularities, but not on the set of ^,x,Vx so that (|2.2.24p and 
(|2.2.7P (resp. (|2.2.1ip and (|2.2.12p ) are satisfied. That set is simply connected, 
so we can modify on its complement and extend it smoothly to f2 x IR x IR^. 
The modification is chosen to depend smoothly on A. 



Require uniform obliqueness: 



(2.2.27) 



Functional independence in upper corners: for i,j = 1,4 and for i,j = 2,3 set 



G 



9p^ 9p^ 



_9p2 9p2_ 

regard it as a function of ^ (including the dependence on Vx(C)) ^nd require 
\\G\\,\\G-^\\<Cd inB5,(ec)nn. (2.2.28) 
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Let T be the se10 of admissible functions so that all of these conditions are 
satisfied. Define !F to be the set of admissible functions such that all of these 
conditions are satisfied with strict inequalities, i.e. replace <,> by <,>, "in- 
creasing" by "strictly increasing" etc. 



[This is the end of Definition l^?^ ] 



The elliptic problem is solved by iteration; is the new iterate, "0 the old one. 
£ defines ifj, as we show later. As always, the iteration is designed so that its 
fixed points solve the problem: 

Remark 2.2.4. liip^tp, then p.2.14p . (|2.2.16p . (|2.2.15p . (|2.2.17|1 and the defi- 
nition of S yield 

(c^/ - Vx^) : V^V- = in n, 

Vx • n = on A and B, 
= X and 
[pVx - p/Vx^) • n = on S, 



L = Vl - e on P 
(we may take closures by regularity (|2.2.ip '). 

Remark 2.2.5. Consider a coordinate system where ^ab = 0. For any point 
on A or B, we can use even reflection of tp across the corresponding boundary 
to obtain a new situation where the point is in the interior. (In or ^b, we 
obtain a new situation with a point at a shock rcsp. quasi-parabolic arc with 
an elliptic region on one side.) The boundary condition Xn = "0" = (due to 
£,AB = 0), for even reflection of tp, implies that tp is across the boundary; 
then necessarily it is also C^'". 

For fixed points ip = ip, standard regularity theory immediately yields that the 
solution is locally analytic (even after reflection). The same technique applied 
to Tp and to solutions ip of linearized equations (here 0, ip and tJj arc reflected) 
yields C^'" regularity. (The same argument applies to S extended by mirror 
reflection across A.) 

Proposition 2.2.6. For sufficiently small e (with bound depending only on Cpt ) 
and rj (depending continuously and only on ip^S.^x)-' 



for all Tp e T , £(V','0') is well-defined for ^j' nearip, and the Frechet derivative 
dC/ dip' {tpi,ip) (of C with respect to its second argument ip' , evaluated at tJ/ — ip) 

-/J 



is a linear isomorphism of C^g" onto Y. 



^The notation does not necessarily imply that is the closure of JF. 
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Proof. The proof is almost identical to |EL07[ Proposition 4.4.6]; the new corner 
between A, B is covered by |Lie88[ Theorem 1.4] in the same way as the other 
ones. □ 

Proposition 2.2.7. 77 can he chosen so that tJj is unique and depends contin- 
uously on ip G J- (both in the C^'" topology) and A. 

Proof. The proof is exactly the same as for |EL07[ Proposition 4.4.7]. □ 

Proposition 2.2.8. For e and rj sufficiently small: for all continuous paths 
t e [0, 1] X{t) in A, Ute(o,i) ('t^J' ^ ^Ht)) °P^'^ Uie[o,i] ({0 x ^Ht)) 
is closcE m [0,1] X C^'"([0,l]2). 

Proof. All conditions on ijj in Definition 12.2.31 arc inequalities which can be 
made scalar by taking a suitable suprcmum or infimum. Then their sides are 
continuous under C|'"([0, 1]^) changes to ip which, by Proposition 12.2.71 means 

continuous in C^'"([0, 1]^) change to (Most inequalities need only C^([0, 1]^).) 

1. Closedness: consider sequences {tn,^n) in Ute[o.i] ({*} ^ •^A(t)) that con- 
verge to a limit (t, tp). 

Let ipn be associated to ipn as in Definition 12.2.31 By continuity (Propo- 
sition , (V'n) converges to a limit ip as well. By continuity of £ in V', 
■0 and A, we have C\(t)['4') '0) = as well. 

Let Sn be defined by as in (j2.2.2p . with s ^ Sn and ip ipn- Then 
by (|2.2.2p . (s„) converges in C^'"[0, 1] as well, to a limit s which satisfies 
((2X2)) itself. 

Most conditions on ?/' are nonstrict inequalities with continuous left- and 
right-hand side, so they are still satisfied by ip- We check the strict in- 
equalities explicitly and in order: 

f^^T^ is imphed by (PX7|) . 

(|2.2.1ip resp. (|2.2.12p resp. (|2.2.13p are implied by (|2.2.10p resp. (|2.2.24p 
resp. (|2.2.9p . by choosing rj sufficiently small. 

All inequalities arc satisfied, so ^ & !F . 

2. Openness: 

same proof, using that all inequalities are strict now, by definition of T , 
hence preserved by sufficiently small perturbations. 

*We make no statement about T being the closure of J^. It certainly contains the closure, 
but it could be bigger, for example if one of the inequalities in Definition 12.2.31 becomes 
nonstrict in the interior without being violated. 
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□ 

Definition 2.2.9. Define JC : T ^ C^'"([0,l]2) to map "0 into as given in 
Definition l2.2.3l but pulled back to (ct, Q coordinates and the [0, 1]^ domain (see 
Definition 12. 2. 3p with the coordinate transform defined by V'- 

2.3 Regularity and compactness 

Proposition 2.3.1. For sufficiently small a E (0,1) and (3 G (1,2), depending 
only on Cd, 5 Lb ■ e, 5o, Cl, ^^x'- 

1. When parametrized in the coordinates of Fiaure \10\. 

\\S\\co,i < CsL (2.3.1) 

and 

||5||p2,. < Cs (2.3.2) 

for CsL = Csl{Cl, 5yx) and Cs ~ C^iCc^ ^vx),' the weight f3 is with respect 
to the endpoints £,At£,c- 

2. For a fixed point ip of K.: 

(a) 112.2. 19]) is strict for sufficiently large C^. 

(b) i2.2.1\) is strict for sufficiently large Cc = Cc(Cd,5Lb-^,CL-,5o,5vtA,5d). 

(c) For K ^n - P ~ {^B, Ub} and all k >0, a' e (0, 1), 

Mc^..'^K)<CcK (2.3.3) 

where Cqk = CcK{d,CL,So,SvtA) is decreasing in d := d{K,PU 
{^ab}) and not dependent on e. 

(d) ip is analytic in Q ~ {S.aBiS.c}! S is analytic except in ^c- 

3. For sufficiently small rj > 0, depending continuously and only on ip, there 
are Sa, > so that for all E T , 

U\y,.^s^m<CK (2.3.4) 
Here, Cic,6a, 5p depend only on Cd, Shb ■ e, (5o, Cl, Svx, 

Proof. The proof is as the one for [EL 07} Proposition 4.5.2], with obvious mod- 
ifications. The only additional problem is the corner in ^ab- This is very easy 
to treat with |EL07( Proposition 5.1.1] because of (|2.2.28|1 and (|2.2.27p for ^is- 
Note that the corner angle in ^ab is bounded away from tt because of the re- 
strictions on ^AB (see Section [2T|) . □ 
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Remark 2.3.2. (|2.3.4p implies in particular that /C is a compact map. ip S 
C^'"([0, 1]^) is mapped continuously into lip € C'^'"^^^" (D,). The latter space 
is compactly embedded in C^'"(0). PuUback to C^'"([0, l]'^) by the aX coor- 
dinates defined by ijj (not ^) may destroy the extra regularity, but preserves 
compactness. 

2.4 Pseudo-Mach number control 

Proposition 2.4.1. For e and S^b sufficiently small, with bounds depending 
only on Sp: if tp €z J- is a fixed point of IC, then i2.2.8]) is strict and 

L^<l-e mQ-P. (2.4.1) 

Proof. 

d(n,p("))>i.e, 

for e small enough. Remember from Definition 12 . 2.21 that b 
fore: 

L^ = l-e<l- WbWco.i ■ d(p(^),p(")) <l-6Lb-b 
e.g. for (5l6 < 1. 

On the shock, we may use p.2.7p combined with |EL07|, Proposition 3.6.1] to 
rule out that L^ + 6Lb-b has a maximum in a point where L < 1 and L > l — S^s, 
with Sls as supplied by loc.cit. Here S^i, has to be chosen so that |(5LhV5| < 5ls 
is satisfied. (Now S^h depends continuously on Sp as well.) 

In addition we can choose 6 Lb so small that SLb ■ b satisfies the preconditions 
of Theorem 1 and Theorem 2 in [EL05aj (where it is called 6) . For Theorem 2 
we use that fe„ = on A and on B. Let Slu be the 6 from those theorems (it 
depends only and continuously on A) . Then + SLh-b cannot have a maximum 
in a point of U A U i? where > 1 — Slu- 

In the corner between A, B, due to regularity the boundary conditions imply 
Vx = 0, so L = 0, so L'^ + SLb ■ b = SLb ■ 6 < 1 for SLb sufficiently small. 

In we use that the shock is pseudo-normal (by the boundary condition Vx ■ 
fiA = which implies Vx-i — for the corresponding shock tangent t since 5*, A 
form a right angle), so Ld = which is uniformly bounded above away from 
1 by a constant depending on Sp, since (j2.2.7p implies uniform shock strength. 

Assume that (|2.2.8p is not strict (or violated). Then L'^ + SLb ■ b has a maximum 
> 1 somewhere. For SLb sufficiently small (no new dependencies) that means 



on P(°). There- 



onp(^\ 
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Figure 15: In this frame P is centered in u/? = 0, i? is horizontal and vj is 
verticaL 

has a maximum > 1 — min{(5Lsi, f^Ls} somewhere. But no matter where in 
12 this occurs, it contradicts one of the cases discussed above. Hence ()2.2.8p is 
strict. 

(|2.4.1[) can be shown in the same manner, by taking 6 = instead, using the 
actual boundary condition L = \/l — e on P and and considering e < (5ls, Slq- 



2.5 Arc control and corner bounds 

The discussion of parabolic arcs is very similar to |EL07[ Sections 4.7 to 4.10]. 
For the convenience of the reader we restate the results using new notation and 
point out some differences in details. 

A new choice of coordinates is convenient (see Figure [T5)) : since self-similar 
potential flow is invariant under translations, we may translate so that vr moves 
to the origin (all other velocities v and coordinates ^ have vr subtracted), then 
rotate clockwise until R is horizontal. In this frame, vi is vertical down and P 
is centered in vr = 0. This means and ^'^e both constant on P, which 
simplifies certain calculations. 

In polar coordinates (r, (j)) with respect to the origin (center of P), P corresponds 
to r = cr ■ \/l — e. 

Proposition 2.5.1. If Cpt < oo is sufficiently large, if Spn > is sufficiently 
small, if e is sufficiently small and Cpy,Cpp sufficiently large, with bounds de- 
pending only on Cpt, then for any fixed point x of IC, h2.2.20\) and \2.2.21\) are 
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strict, and 

\P~ Pii\<Cppe^/^ and (2.5.1) 
\v- vr\ < Cp^e^l'^ on P. (2.5.2) 

Proof. The proof is as for [EL07[ Proposition 4.8.1], with obvious modifications. 

□ 

If ()2.2.4|) is satisfied, but not in its strict version, then 77^ = 77^ or 77^ = rj^ 
(where are as defined in Definition 12.2.31 after (|2.2.4p ). Each of these two 
cases must be ruled out. 

Proposition 2.5.2. For e sufficiently small: for any fixed point ijj € of JC, 
the lower bound in \2.2.J^ is strict: 

Vc > Vc 

Proof Same as for |EL07[ Proposition 4.10.1]. □ 

Proposition 2.5.3. Consider rjc = jy^. For sufficiently small e, there is an 
a > so that 

1. tjj + does not have a local minimum (with respect to CI) at P U {^b}, 
and 

2. a shock through with upstream data vj and pj and tangent (1, -^^) has 

yy > 0. 

Proof This foUows as in Propositions 4.10.2, 4.10.3 and 4.10.5 of [EL07| . □ 
Only the final upper bound requires some adaptation: 

Proposition 2.5.4. Let x ^ ^ be a fixed point of K.. For Cn sufficiently large 
and for e > sufficiently small, the upper part of {2.2.4-\ l is strict: 

VC < vi- 

Proof. Again, consider the coordinates of Figure [TSl 

By Proposition 12.5.31 V' + o,^ cannot have a local minimum at P U {£,b}- For 
Vc = »7c, we have {ip + a^)^ = -0,, > in by [EL07| (4.9.8)] (for sufficiently 
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small e), so the minimum cannot be in either (note that the domain locally 
contains the ray downward from the corner). 

On the shock (excluding endpoints): let ^ i— > s(^) be a local parametrization of 
the shock . + a£, = + a^, so 



dt{^ + aO = dtii'' + ai) = VI ■ t + 



VjS^ + a 



(1 + 4)1/2 (l + s|)l/2- 

For a local minimum at the shock we need dt{ip + a^) = 0, so 

a 

= y- 

A global minimum, in particular < '0(Cc)+a^c, additionally requires that (as 
well as the rest of the shock) is on or below the tangent through the minimum 
point, because ij}^ and thus -0^+0^ are decreasing in -q. By Proposition l2.5.3l the 
shock through with that tangent has > for rjc ~ Vc- ^''^ minimum 
point the tangent has same slope but is at least as high, so the shock speed is 
at least as high, so = i/j^ = [t(j + a^),, there is at least as high, in particular 
> too. But that contradicts a minimum (the ray vertically downwards from 
any shock point is locally contained in fl, by (j2.2.26p ). Hence ip + cannot 
have a global minimum at the shock. 

The equation (2.2.5) yields 

(a^ is linear), so the classical strong maximum principle rules out a minimum 
in the interior (unless "0 + is constant, which means we are looking at the 
unperturbed solution which has rjc = i]q < Vc)- 

On B, the boundary condition V'n = Xn = implies {ip + a^)„ = a^„ > (the 
slope of B in the frame of Figure [T5|) is always nonnegative), so the Hopf lemma 
rules out a minimum of "0 + at B. 

On A the boundary condition Xn = yields ipn = ^-n — (,ab -^a > (see Figure 
[T5|) . This is actually i/jn > 0, except in the special case where (in the notation 
of Definition I2.1.2p (^ab = "^/j which is allowed only if 7 = 1 and i]^ — rj^: 
the "unperturbed" case. In that case, the proof of Proposition 12.8.11 shows that 
only the unperturbed solution (Figure [T5)) can solve the problem. Its corner is 
exactly in the expected location, so that Vc = Vc ^ Vc ■ '-' 



2.6 Velocity and shock normal control 



Proposition 2.6.1. If C'i,tF(.TCynA o,re sufficiently large (bounds depending only 
on Cpt), if Csn is sufficiently large (bound depending only on CytRiCynA), if 
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VI 




Figure 16: A maximum of requires negative curvature, causing a contradic- 
tion 




Figure 17: Left: mirror-reflect Figure [T^] across A and rotate around the origin. 
Right: setting of Figure [HI left. 



February 5, 2008 



39 



2.6 Velocity and shock normal controB CONSTRUCTION OF THE FLOW 



e is sufficiently small (bound depending only on Csn), and if 5cc is sufficiently 
small, then for any fixed point ip € T of IC, the inequalities I12.2.23\) . i2.2.22\) . 
\2.2.26\) and i2.2.5\) are strict. Moreover 

\Xt\>S^t onSnBs.iCc), (2.6.1) 
for some constants 5^t, > 0. 

Proof. 1. For (|2.2.23p : consider the coordinates of Figure [T51 where tR = 
(1,0). Let ^ ^ s(^) parametrize S (the shock normal bounds (|2.2.26p 
show that S is nowhere vertical in these coordinates,, for sufficiently small 
e, bound depending on Csn)- Assume that v -tif = attains a positive 
global maximum (with respect to VL) in a point at S (i.e. on the down- 
stream side). Since u/ = (0, Wj) with < 0, this means < in 
(because < 0), i.e. 54(^0) < (see Figure [TG]). 

sj(^o) can be expressed as a continuous function of w^(^o) and ^o- The set 
of possible ^0 is contained in the set of possible shock locations which is 
prc-compact. Therefore if = C^tR ■ e^^^ hi ^0 £ S, then 

.s?(eo) < • e'/' (2.6.2) 

where Csi = Csi{CutFi.) > is uniformly increasing in CytA- 

For a constant-state solution (|2.2.23p is immediate. Otherwise, since S 
and V' are analytic (Proposition 12 . 3 . I|) . we can apply [EL07[ Proposition 
3.5.1] with w = (1,0), which yields that curvature k < 0, i.e. sjj > 0, 
in ^0- Therefore s^{£_) < s^(^o) for C < Co near ^o- On the other hand, 
> in C^ since the boundary condition Xn — requires the shock to be 
perpendicular to the wall A; in particular s^{^a) > > 55(^0) by (|2.6.2p . 
(In this choice of coordinates, A is either vertical or has negative slope, 
since we require it to form right or sharp angles with R, by choice of ^ab 
in Section O) 

Therefore wc can pick S,a € iS,A, Co) maximal so that sj(Ca) = sc(Co)- Then 
sdO < s?(Co) for C € (Ca,Co), so by integration 

5(Ca) >5(eo) + sdCo)-(Ca-eo). 

But that means the shock tangent in Ca is parallel to the one in Co but 
higher, so the shock speed a := £^ ■ n is smaller. By jEL07[ (2.4.19)], that 
means is smaller, whereas is the same (parallel tangents), < 0, so 
is bigger. Contradiction — we assumed that we have a global maximum 
of v'^ in Co- 

[EL07[ Propositions 3.3.1 and 3.4.1] rule out local maxima of in fl and 
at B, where we use that x is analytic and that (1,0) is not vertical, i.e. 
not normal to B. 
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At A: if A is vertical, then the boundary condition requires = £,a < 0; 
if A is not vertical, then (1,0) is not normal, so [EL07[ Proposition 3.3.1] 
applies again. 

In ^AB, the two boundary conditions combine to yield v = ^ab, so = 
Ub < 0. 

In ^A, > (see above) yields < 0. 

On P we can use (|2.5.2|) with vf^ — 0, increasing CytR to > Cp^ if necessary 
(this makes CytB, depend on Cpt as well). 

All parts of H, are covered; p.2.23p is strict. 

2. For (|2.2.22p : consider the coordinates of Figure [T7] left. There, = (1, 0), 
so we need to show = v-ua < CmtA'^-^^'^ ■ On A, the boundary condition 
yields = 0. B is never vertical, so |EL07[ Proposition 3.4.1] rules out 
extrema of at B. [EL07[ Proposition 3.3.1] does not allow cxtrcma in 
n. At P, (PX^ yields ii^ = + 0{e^/'^)] note that < in these 
coordinates. At 5, we can use the same curvature argument as for v ■ tp, 
except that we now use > in rather than ^a- Altogether we obtain 
a contradiction again, if CynA is sufficiently large, depending only and 
continuously on Cpt- 

3. Consider the coordinates of Figure [TSl The slope of some shock passing 
through a point ^ is uniquely determined by (and continuous in) ^ and 
v^, with sgns^ = — sgnw"^ (since vi = (0,z;f), < 0). The set of possible 
shock locations ^ is pre-compact, so p.2.23|) implies 

sup Z(n, Hp) < Csn ■ e^^^ 

where Csn = Csn{CytR)- 

Analogously we argue that (|2.2.25[) implies 

supZ(Ci,n) < Csn ■ e^^^ 

where Csn = Csn(Ct)ti?,, Cim^) now. (|2.2.26p is strict with these choices. 

4. These shock normal bounds also imply (|2.2.5[) is strict, for 5cc > and 
e > sufficiently small(er), with e bound depending only on C'sn- 

5. Near each corner the shock normal bound bounds n away from the ^ 
direction, so \xi\ > Sy^t and therefore (|2.6.ip for some S^^t. 

□ 

Proposition 2.6.2. 1. If 6sb is sufficiently small, then \2.2.3\] is strict. 
2. There is a constant 6ps > so that 

Pd> pi + Sps at S (2.6.3) 



February 5, 2008 



41 



2.6 Velocity and shock normal controB CONSTRUCTION OF THE FLOW 



Proof. 1. Consider the envelope E defined in Section [??T] The parameter set 
A (see Definition 12. 1.2p has been chosen so that for any X G A, E passes 
from ^'"^ to A without meeting B or the circle (with radius c/ centered in 
vi). Since A has also been chosen compact, E is in fact uniformly bounded 
away from B and the circle. 

E starts in ^p"''; let E' be the counterclockwise envelope (Definition [TT73]) 
starting in instead. E,E' are solutions of an ODE (|1.7.ip . so they 
depend continuously on the initial point. Hence for sufficiently close to 
Cc^'' ^ i-C- by (|2.2.4p for sufficiently small e (with upper bound depending 
only on the choice of A). E' is also uniformly bounded away from B and 
the circle. 

Now we can apply the argument displayed in Figure [S] right: |^ — w/j is 
r in the polar coordinates used in Section 11.71 Let E' and the shock S 
be parametrized hy (j> i-^ rs[4>) resp. 4> ^ i^E'i4>)j with (j) G [(pCi'pA], 4>c 
corresponding to the ray from vj through S^c a-nd <j)A to the ray from Vi 
containing A. rs{4>c) = TE'{(t>c) because S and E' both pass through £^c- 
By (12211), id < 1 at S. Therefore, Proposition[L72]yiclds rsi^) > rE'lffy 
for all (j> > (pc- Hence topologically S is separated from B and the circle 
by E' , so it also has uniformly lower bounded distance from them. In 
particular (|2.2.3p is strict, for sufficiently small Ssb (depending only on 
the choice of A, but not on any other constant). 

2. If S vanishes in some point ^, then Ld = — \£, — vj\/ci which — since S 
has uniform distance from the circle — is uniformly bounded below away 
from 1. However, this contradicts (|2.2.9p . The shock cannot vanish; on the 
contrary, by continuity the shock has uniformly lower-bounded strength. 
That implies (|2.6.3p . for sufficiently small 6ps- (Again, it depends only on 
A, not on the choice of other constants.) 

□ 

Proposition 2.6.3. If Sp and e are sufficiently small (with bounds depending 
only on Cpt), then for any fixed point ip £ T of JC, the inequality 12.2.T\) is 
strict. 

Proof. By Proposition 12.3.11 V and hence s are analytic. Thus we may use 
|EL07[ Proposition 3.2.1] which rules out minima of p in f7 and (using Remark 
[2X5)) at A or B. 

Consider the coordinates of Figure [121 In ■Ca, the first shock condition is 
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(|2.2.25p implies 

So in S^AB — 0, since Vx = by boundary conditions on A, B and regularity; 

P = ^"'(-X - \\^X?) = ^"'(-^) = PI + 5pAB 

for some constant 5pAB > depending only on the parameters A; note that tt 
is a strictly increasing function for any 7 > 1. We can pick 5p < 6pAB so that 
p < Pi + Sp is not possible in ^ab ■ 

On P we know p up to a small constant, by ()2.5.ip . so we can choose 6p even 
smaller so that p < pi + Sp is not possible at P. 

By (|2.6.3p . p at is uniformly bounded below away from pj. Hence, for Sp 
sufficiently small, p cannot have a global minimum close to pi at S. 

We see that for sufficiently small Sp and e, depending continuously on Cpt (and 
A), ([QJ]) is strict □ 

Proposition 2.6.4. IfS^tA, SynB ande are sufficiently small (SytA, SvnB bounds 
depending only on Sp,Csn, e hound depending only on Cpt), and if Scc is suffi- 
ciently small, then for any fixed point ip (z of tC, the inequalities 12. 2. 25\) and 
\2.2.24\j are strict. 

Proof. Consider the coordinates of Figure [T7l right, where v ■ ub = —v^. (|2.2.7p 
implies that the shock is uniformly strong. By (|2.2.26p . the shock normal n is 
everywhere downwards and uniformly not horizontal. Thus > Vj + SynB at 
S for sufficiently small SynB, depending only on Sp and Csn- 

|EL07[ Proposition 3.3.1] rules out local maxima of in fl. 

If has a local maximum at A, then A must be horizontal ( |EL07| Proposition 
3.4.1]), but by construction it is not. 

On B the boundary condition implies = x„ = X2: so = ip2 = Vab = 0. 

At P we can use ()2.5.2|) with Vp = 0> v\ to obtain > Vj ii e is small enough 
(depending on Cpt). 

Altogether we have that (|2.2.24p is strict if SynB is small enough. 

The arguments for (|2.2.25[) are analogous, looking at Figure [TTl left coordinates 
instead: the shock S is nowhere vertical (by ()2.2.26p ). so > + S^tA at S 
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for sufficiently small SytA- If B is not horizontal, then the direction (0, 1) is not 
perpendicular to it, so |EL07[ Proposition 3.4.1] rules out a local extremum 
at B; if B is horizontal, then = Xn = X2, so ~ ■02 = Vab = on it. A 
is always vertical, i.e. never perpendicular to (0,1), so by [EL07[ Proposition 
3.4.1] no extremum is possible at it. In ^ab = 0, the boundary conditions 
combine to -u = 0, so = > + S^tA if S^tA is small enough. At P we can 
use (|2.5.2p again to obtain > v^ — Cpt ■ e^^^ > Vj (using > and for 
e sufficiently small, with bound depending only on Cpt). [EL07[ Proposition 
3.3.1] rules out interior extrema of v'^ . Hence (|2.2.25p is strict if SytA is small 
enough. □ 



2.7 Fixed points 



Proposition 2.7.1. For So sufficiently small, with bounds depending only on 
Sp and Cl, for Cd resp. Sd sufficiently large resp. small, with bounds depending 
only on Sp and Cl , and for e sufficiently small, with bounds depending only on 
Cpt, Cl and Sp.- 

If X is a fixed point of JC, then {2.2.2T\l and 12.2.28]) are strict. 

Proof. Compared to |EL07[ Proposition 4.13.1], the only new case is a corner 
between two walls, A and B. The corner angle is bounded away from and tt by 
constants depending only on the parameters A. (Note that S,ab in Section 12.11 
has been lower-bounded uniformly by in Definition 12. 1.21 so that A, B are 
uniformly not parallel.) on A and B is their respective normal, so (|2.2.28p is 
obvious. □ 

Proposition 2.7.2. // the constants in \2.1.2]) in Definition \2.2.S\ are chosen 
sufficiently small resp. large: for any A G A, /Ca cannot have fixed points on 

Proof. Let x G be a fixed point of K.. We show that every inequality in the 
definition of T is strict, so x G J^. 

(|2.2.ip and (|2.2.19p arc strict by Proposition HXH 

(|2.2.3p is strict by Proposition [2X1 

(|2.2.7p is strict by Proposition [Mill 

A fixed point satisfies ?/' = V': so ||^/' — i/'H = r/('0) > cannot be true. ()2.2.10p 
is strict. 
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(|2.2.8p strict is provided by Proposition 12.4.11 

Due to Proposition 12.4.11 = 1 — e on eaeh point of P, so we are in the 
situation of Section [2T5l and [EL07| Section 4.7 etc]. Proposition 12.5.1) shows 
that (|2.2.20|1 and ()2.2.2ip are strict. 

(|2.2.22p and (|2.2.23p are strict by Proposition [2XT1 

(|2.2.24p and (j2.2.25p are strict by Proposition [Mil 

Propositions 12.5.^ and [^5.41 rule out ijc = Vc ^ S^^e if S is small enough, so 
((2:2:4)) is strict. 

(|2.2.5p is strict by Proposition EXH 

(|2.2.ip yields a trivial upper bound on the density in Q, hence downstream at 
the shock. 

(|2.2.26p is strict by Proposition [2:0 

Proposition 12.7. II shows that (|2.2.27p and p.2.28p are strict. 

All inequalities are strict, so ip G □ 



2.8 Existence of fixed points 

We determine the Leray-Schauder degree of IC on J- for a particular choice of pa- 
rameters A: the unperturbed problem (see Figure [15)1 . featuring a straight shock 
separating two constant-state regions = 77^, ^ab = fff in the coordinates 
of Definition [2T2I), for 7 = 1. 

Proposition 2.8.1. For sufficiently small e: 

For 7=1, rjQ — rj*Q and £^ab — vff, K, has nonzero Leray-Schauder degree. 

Proof. We can use reflection across A (Remark 12. 2. 5p to obtain the problem of 
Propositions 4.14.1 and 4.14.3 in |EL07| . The resulting iteration JC is almost 
the same as in loc.cit., except for minor differences in the coordinate transform 
from {cr,C,) £ [Oj 1]^ (fixed domain) to ^ coordinates (see Definition 12.2.31 as 
compared to |EL07[ Definition 4.4.3]). The proofs of |EL07[ Propositions 4.14.1 
and 4.14.3] carry over without any change to show that the present problem has 
nonzero Leray-Schauder degree. □ 
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1 
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^■^ constant 
/■ density, 
velocity 

\ 
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Take 9 T 90° 
reflection point "f 


oo 



Figure 18: The unperturbed case: a straight vertical shock R. In this case there 
is no reflection point and no incident shock. 

Proposition 2.8.2. For sufficiently small resp. large constants in K. 
has a fixed point for all X Cz A. 

Proof. The proof is identical to the one of |EL071 Proposition 4.15.1], except for 
the definition of A (Dcfinition l2.1.2|) : we use the known Leray-Schauder degree 
in (7, ric, £,ab) = (1, Vc^ ^'b) from Proposition |2A1J □ 



2.9 Construction of the entire flow 

Proof of Theorem[l[ For all pi,cj,Mi G (0, cxd) and for each choice (in Def- 
inition I2.1.2|) of 7, 77^ and we obtain a separate parameter set A. For 
sufficiently small constants in ()2T.2p . Proposition 12 . 8 . 2] vields fixed points tp for 
all A e A. Note that there is no lower bound on e, except that a, (3 etc. may 
change as e 4 0. 

By Definitions^ Remark[M!l Proposition l^XT] and ((03)) . the fixed points 
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I 
I 




Figure 19: The expected and actual parabolic arc (P*^'^) and P^*^^) differ by 
curve of length 0{e^/^) (by ((2:2:41) ) 



February 5, 2008 



41 



2.9 Construction of the entire Bow 2 CONSTRUCTION OF THE FLOW 



satisfy 





= 


in n^'\ 


(2.9.1) 






and 


(2.9.2) 


\p- pr\ 


= 0(ei/2) 


and 


(2.9.3) 


|VV - Vr\ 


= 0(ei/2) 


on P(^\ 


(2.9.4) 


X 




and 


(2.9.5) 


(pVx- Pi^X^) -n 


= 


on 5*, 


(2.9.6) 


Vx • n 


= 


on A(JB, 


(2.9.7) 




^ 0(ei/2) 




(2.9.8) 



where the O constants are independent of e. For regularity, Proposition I2.3.11 
yields 

\mco^0'^J < Ci, (2.9.9) 
IIV^IIc^,= (W^')'l^lc-=(W=')^^2(rf) (2.9.10) 
where d := d{K, P'^'^ U {Ub}) > 0. 
for constants Ci and C2{d) independent of e. 

Now consider those parameter vectors A that arise from the situtation in The- 
orem [1] i.e. so that there is an incident shock Q meeting i? in a local regular 
reflection. We extend ip from above to a function -0^'^ defined on all of V as 
shown in Figure [T9l set p = p^, v ^ vji in the region enclosed by R shock, B 
and set p = pq, v = vb.,q in the region right of the Q shock and p = pi, 

V = VI in the remaining area. In each of the four regions, tp^'^'^ is a strong solu- 
tion of self-similar potential flow, so we can multiply the divergence-form PDF 
|EL07[ (2.2.3)] with any test function -& € C^(V) and integrate over all region 
to obtain a sum of boundary integrals of the type 

pVx • n ds 

M 

where M are various curves; Vx and p are limits on one of the sides of M . 

The symmetric difference of P^"'^ and P^"^ has length 0{e^/'^) (by (Hm, so 
since Vijj and are bounded in each region (uniformly in e, by (|2.9.9p ). the 
boundary integral over the difference contributes only 0{e^/^). The difference 
of the integrals on each side of P*(^'nP(") are ©(e^/^) due to ^^T^ and (PTO)) . 
The integrals over A, B vanish due to (|2.9.7p . Finally, the integrals on each side 
of S'(<^) cancel due to ((2X5)) and ^MM- Ahogether: 

pWVx^') • Vd - 2p^'^ d^= 0(6^/2). (2.9.11) 

V' 
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C^-°' with fc + a > 1 is compactly embedded in C°'^, so by (|2.9.10p with a 

diagonalization argument, for every compact K d V — {£,ab} — P we can 
find a sequence (e^) J, so that ip^'^'''^ converges to 

^(0) C°'^{K). Moreover 
p*^"^^ and Vx^^^ are bounded on V imiformly in e, so we may take e | in (|2.9.1ip 
to obtain 

/ p(o)Vx^"^ • - 2p('')?9 0. (2.9.12) 
Jv 

In addition, (|2.9.5|) and (|2.9.2|) combined with (|2.9.9D show that 

V'^"^ e C{V) (2.9.13) 

Finally, by construction of 

p(")(.a,«^°H^O-^''"i' C^l^'' as.^oo, (2.9.14) 

i.e. their limits on rays to infinity are exactly as for the initial data in Figured 
This means the limit approaches the initial data as t I 0. 

(|2.9.12p . (ElJl, (|2.9.13P and (|2.9.14p show that 0(i,f) := ip'-°^t-^x) defines a 
solution of 11133), ((OSl), (fLO)) and ([03)1 . 

By taking 7 f oo, 77^ j and i C_bb, we obtain a solution for euery 

7 G [1, 00), G [77^, 0) and ^ab G {£,eB, v%}. (in the cases 7 > 1 and 77^ = 77^, 
we may use that rj^ approaches t/^ as e | 0). 

As mentioned (Remark l2.1.ip . this exhausts all cases covered by the conditions 
of Theorem [TJ The proof is therefore complete. □ 

Remark 2.9.1. In addition to mere existence we obtain some structural infor- 
mation in the proof: 

1. The solution has the structure shown in Figure [4] left, with pscudo-Mach 
number L > 1 in the /, R, Q regions, L < 1 in the elliptic region Q. 

2. The solution has constant density and velocity in each of the /, R, Q re- 
gions. 

— *(o) . 

3. The solution is analytic everywhere except perhaps at P and in £,ab 
and. of course, the shocks. 

— *(o) 

4. The curved shock is analytic away from A and P and Lipschitz overall. 

5. Density and velocity are bounded. 
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It is expected that density and velocity are at least continuous. However, the 
methods developed in |EL07j yield boundedness everywhere, but continuity only 
away from P . Note that P can not be a classical shock with smooth data on 
each side, because the one-sided limit of L on the hyperbolic side R of P* is 
= 1 everywhere (> 1 is needed for positive shock strength). 

Some additional structural information: 

1. The possible (downstream) normals of the curved shock are between n^j 
and tA (counterclockwise). 

2. The shocks are admissible and do not vanish anywhere. 

3. In the elliptic region, < vf and > (in Figure 2] left coordinates). 

4. In the elliptic region, the density p is greater than pj. 

Additional information can be obtained from the inequalities in Definition 1 2. 2. 31 
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